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Abstract 



We investigate the prospects of indirect searches for supersymmetric neu- 
tralino dark matter. Relic neutralinos gravitationally accumulate in the Sun 
and their annihilations produce high energy neutrinos. Muon neutrinos of this 
origin can be seen in large detectors like AMANDA, IceCube and ANTARES. 
We evaluate the relic density and the detection rate in several models — the 
minimal supersymmetric model, minimal supergravity, and supergravity with 
non-universal Higgs boson masses at the grand unification scale. We make re- 
alistic estimates for the indirect detection rates including effects of the muon 
detection threshold, quark hadronization, and solar absorption. We find good 
prospects for detection of neutralinos with mass above 200 GeV. 
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I. INTRODUCTION 



Average matter and energy densities in the Universe are commonly described in 
terms of density parameters (fij = Pi/p c ), where 

p c = 3H 2 /(8nG N ) ~ 1.88 x 10~ 29 /i 2 g/cm 3 , (1) 

is the critical density to close the Universe, h = H /(100 km sec -1 Mpc -1 ), H is the Hubble 
constant, and is Newton's gravitational constant. Recent measurements of the Hubble 
constant are converging to h ~ 0.6 — 0.7 

Studies of clusters of galaxies and large scale structure indicate that the matter density 
should be at least Qm ^ 0.2 ||. The matter density inferred from observations of cluster X- 
ray ||, supernovae ||, and the cosmic microwave background radiation (CMB) anisotropy 
data [0-0 is Qm ~ 0.3 — 0.4. The baryon density inferred from Big-Bang Nucleosynthesis is 
Q b h 2 ~ 0.020 ± 0.002 |K| which is in good agreement with the value obtained from recent 
analyses of the CMB power spectrum The luminous matter in the Universe has a very 
small relative density fl^ ^ 0.01 ||11||. Therefore, the evidence for substantial dark matter is 



compelling. Moreover, the dominant component must be cold dark matter to be consistent 
with large scale structure |12[]. The most attractive cold dark matter candidates are stable 
weakly interacting mass particles (WIMPs) produced in the early Universe. In our analysis, 
we take Q x = Qm — to be the WIMP density and conservatively consider 

0.05 < n x h 2 < 0.3. (2) 

as the cosmologically interesting region for the cold dark matter density . 

In the early Universe, the temperature (T) was much higher than the WIMP mass {m x ) 
and WIMPs existed abundantly in thermal equilibrium with WIMP annihilation into lighter 
particles balanced by pair production. When the Universe cooled and the temperature fell 
below m x , the WIMP density became suppressed by the Boltzmann factor exp(— m x /T). 
As a result of the expansion of the Universe, WIMPs dropped out of equilibrium with other 
particles at some point, leaving a relic abundance. 

From simple dimensional analysis, the WIMP annihilation cross section is about 

a ~ a 2 /m x , (3) 

and their relative density is 

n x h 2 = m x n x /p c ~ (3 x 10~ 27 cm 3 /s)/ < av >~ (3 x 10~ 27 cm 3 /s)(m 2 /a 2 ), (4) 

where a is the fine-structure constant. For a WIMP mass of the order of the weak boson 
mass, m x ~ mw, the WIMP density is O x ~ 1. Thus, a cosmologically interesting fl x is 
natural in theories with a weak scale WIMP mass. 

In super symmetric (SUSY) theories with a conserved _R-parity[], the lightest SUSY parti- 
cle (LSP) is stable and is thus an attractive WIMP candidate |T1-|TT5|. Most commonly, the 



R = +1 for particles in the Standard Model as well as Higgs bosons and R = — 1 for their 
superpartners. 
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LSP is the lightest neutralino which is a superposition of the supersymmetric partners of the 
photon, the Z boson and the neutral Higgs bosons. The gaugino and higgsino composition 
of the neutralino is an important factor in the neutralino annihilation cross section. If the 
neutralinos are more gaugino-like, they annihilate mostly into bb and rf; if the neutralinos 
are higgsino-like, most of them annihilate into gauge bosons and ti. 

WIMPs in the halo of our galaxy lose energy when they pass through the Sun and 
scatter from nuclei. When their velocities fall below the escape velocity from the Sun, 
they become gravitationally trapped and accumulate. Annihilations of these WIMPs are a 
potential source of high energy neutrinos that can be detected via Cerenkov light by large 
under-ice or under-water photomultiplier arrays in neutrino telescopes, such as AMANDA 



17 , IceCube 18 and ANTARES [19 



In this paper we investigate the prospects of observing muon neutrinos produced by 
WIMP annihilation in the Sun. Due to cosmic ray backgrounds, only the signals from the 
muon neutrinos at or below the horizontal detection can be detected. In Section II, we 
estimate the general expectations for neutrino signals from WIMP annihilation, the effect 
of the muon detection energy threshold, and the backgrounds to the signal. In Sections 
III, IV and V, we evaluate the indirect detection rate within the framework of the minimal 
supersymmetric model (MSSM), the minimal supergravity model (mSUGRA), and super- 
gravity unified models with non-universal Higgs boson masses at the grand unified scale, 
respectively. In Section VI and VII, we summarize the prospects of indirect search and 
our conclusions about indirect detection of neutralino dark matter in large ice or water 
experiments. 



II. INDIRECT DETECTION OF WIMPS 

When our galaxy formed, the cold dark matter WIMPs were clustered with the luminous 
matter and formed a significant fraction of the galactic matter density [pOj , pl 



Px ~ 0.3 - 0.5 GeV/cm 3 (5) 
as deduced from rotation curves. The distribution of WIMPs in the galactic halo is usually 



approximated as an isothermal sphere with an average velocity p2|-p4|j 



v x ~ 220 - 300 km/sec. (6) 

The density of WIMPs that are gravitationally trapped in the Sun builds up until equilibrium 
is reached when their annihilation rate is half of the capture rate (the factor of 2 is due to the 
Majorana nature of WIMPs). Since the WIMP mass is naturally in the 100 GeV-1000 GeV 
mass range, they annihilate mostly into weak bosons and heavy quarks, which subsequently 
decay into final states with high-energy neutrinos [^,^]. The high-energy muon neutrinos 
can be detected in large ice or water detectors by the Cherenkov light from the muons. The 
muons in these events are aligned along a direction pointing back to the Sun. The angular 



spread of the muons from this direction is 6^ ~ 1.2°/y£' At (TeV). The atmospheric neutrino 
events in a pixel containing this signal are the backgrounds. The measurement of muon 
energy might be used to (i) focus the search window around the direction from the Sun 
in order to reduce the background and (ii) estimate the neutrino energy from the angular 
spread. 
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A. Indirect Detection Rate 



In this subsection, we make a rough estimate of the indirect detection rate for a generic 
WIMP dark matter under the following assumptions^ f25fl : 

(i) the measured galactic halo density is mostly associated with WIMPs and their flux is 

4 cm 3 

X = n x v x ^ ^^-^ x (3 x 10 cm/s) 

1-2 x 10 7 _ 2 _ x 

-cm s , (7) 



m x (GeV) 

where m x is the WIMP mass; 

(ii) the WIMP-nucleon cross section is approximately given by the dimensional analysis 
result 

a( X N) = a DA ~ (G F m 2 N ) 2 -^- = 6 x 1(T 42 cm 2 ; (8) 

m w 

(iii) the WIMPs annihilate about 10% of the time into neutrinos via the channels XX ~ > 
W + W~ or QQ where Q is a heavy quark (bottom or top). 

The further physical quantities needed to calculate the indirect detection rate for neu- 
trinos from WIMP annihilation in the Sun are as follows: 

(i) the capture cross section in the Sun (a ), given by the product of the number of 
target nucleons in the Sun and the elastic \N cross section (uei) and a focusing factor 
(/ ~ 10) that is given by the ratio of kinetic and potential energy of the WIMP near 
the Sun @. 

a Q = f [1.2 x 10 57 ] a n . (9) 

(ii) the neutrino flux at the Earth from WIMP annihilation in the Sun 

<f> v = B x <f> x a Q /4ird 2 , (10) 

where d = 1A.U. = 1.5 x 10 13 cm and B x is the branching fraction to neutrinos 
(~ 10%). Thus we find 

3 x 10~ 5 _ 2 _ x 

cm s 1 . (11) 



m x (GeV) 



In this section we take p x = 0.4 GeV/cm 3 and v x = 300 km/s. 



4 



(iii) the probability of neutrino detection [26], given by 

P = N A a u R^ ~ 2 x 10~ 13 [m x (GeV)] 2 , (12) 

where Na — 6 x 10 23 is the Avogadro number, o v = 0.5 x lCT 38 ^ (GeV) cm 2 is the 
average of the vN and vN cross sections, and ~ 500 cm x E M (GeV) is the muon 
range. For the kinematics of the annihilation and decay chain 

xx - w + w~ 

the average neutrino energy is about < E v >~ |^ x . For the decay chain into heavy 
quarks 

xx - QQ 

the average neutrino energy is about < E v >~ jm x . 

(iv) the number of neutrino events per year per unit area (m 2 ) detectable in a neutrino 
telescope from \X annihilation in the Sun 

dN m /dA = (f) u P~2x 10" 6 m x (GeV) (year)" 1 (m 2 )" 1 . (13) 

For an effective area of 10 4 m 2 , the indirect detection rate is approximately 

dN m /dA - 2x 10" 2 m x (GeV) (10 4 m 2 )- 1 (year)- 1 . (14) 



B. Neutrino-generated Muons 

Neutrino telescopes detect muons from the charged-current reaction 

^ + N -> /x + X, (15) 

where N is an average nucleon in the Earth below the detector. The effective detection 
volume is the product of the detector area and the muon range in rock (R^) which increases 
with energy. By observing muons produced by charged current of v^ipp) interactions in 
the rock below the detector, the effective volume is enhanced. The characteristic range 
of TeV muons is about one kilometer in rock, so large effective volumes are realized at 
TeV energies. This technique applies only for upward going muons entering the detector 
from below the horizontal, because of the large backgrounds from atmospheric muons that 
dominate neutrino signals from above the horizontal. 
The average muon energy loss rate is 

(jjf) = - (3(E) x E, (16) 
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where X is the thickness of material in g/cm 2 . The first term represents ionization losses, 
which are approximately continuous, with a ~ 2 MeV g _1 cm 2 . The second term includes 
the catastrophic processes of bremsstrahlung, pair production and nuclear interactions, in 
which fluctuations play an essential role. Here (3 ~ 4 x 10~ 6 g _1 cm 2 . The critical energy 
above which the radiative processes dominate is 

E cr = a/P « 500 GeV. (17) 

To treat muon propagation properly when > E cr requires a Monte Carlo calculation of 
the probability P surv that a muon of energy E^ survives with energy > E™ m after propagating 
a distance X ||27|| . The probability that a neutrino of energy E u on a trajectory through a 
detector produces a muon above the threshold energy at the detector is P7|-|2T3[| 

P v (K,E™ n ) = N A J o dE^iE^E^R^E^E™), (18) 

where 

R cS = / dXP surv (P M ,Pf n ,X). (19) 

J 

The flux of ^-generated muons at the detector is given by a convolution of the neutrino 
spectrum (j) u with the muon production probability ( |I~8"D as 

MK' m ' 9 ) = I dE u P u (E u ,E™ n )expi-cr tot (K)N A X(6)}ME u ,6). (20) 

M 

The exponential factor in the integrand accounts for absorption of neutrinos along the chord 
through the Earth, X(6). Absorption becomes important for a(E u ) > 10~ 33 cm 2 or E v > 
10 7 GeV. The event rate is calculated by multiplying Eq. ( ^0|) with the effective area of the 
detector. 

The probability P U (E U , 0) can be approximated with a power law 

P_^1.3xl0- 6 (^) 2 - 2 (21) 

for E v between 1 GeV and 1 TeV. This energy dependence reflects the neutrino cross section 
[a v oc E) and the effective muon range (R^ oc E) in Eq. (0). Above 1 TeV, the effect of 
the W propagator slows the growth of a v and the muon range also makes a transition from 
linearity to a more gradual energy dependence. 



C. Atmospheric Neutrino Backgrounds 

The background to indirect detection (ID) is determined by the flux of atmospheric 
neutrinos in a pixel region of the detector that subtends the Sun. The number of upward 
atmospheric background events per year in a 10 4 m 2 detector is ~ 10 2 /P M (TeV). The size of 

the pixel region is determined by the angle between muon and neutrino (~ 1.2° / ^/P^(TeV)). 
Using the kinematics P M ~ m x/^ f° r the WIMP signal from annihilation to PU-pairs, we 
obtain 
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B 



ID 



10 2 /^(TeV) 



2tt, 



1 2° 7r „ 



1.1 x 10 5 

m\ (GeVy 



10 4 m 2 1 year 1 . 



This expression is only valid for for ~ ttl x /4 > 100 GeV. More generally, the backgrounds 
are given in Table I. 

TABLE I. Number of background events due to atmospheric neutrinos. 



£ M (GeV) 


N(Background) 
in 10 4 m 2 


Number of pixels 
of solar size 
in 2-ir 


N(Background) 
per 10 4 m 2 
per pixel, per year 


10 


3200 


140 


23 


100 


1060 


1.4 x 10 3 


0.8 


1000 


110 


1.4 x 10 4 


8 x 10~ 3 



For large m x the signal to background ratio is 

(—) = f ro * 7.2 x 1Q-X (GeV) a, 
KBJm dB m /dA x V ; 

Table II shows the expected signal and signal to background estimates for several values of 
m x . Thus high-energy muons pointing at the Sun are a viable signal for sufficiently high 

m-y. 



TABLE II. Event rates and signal to background (Ns/Nb) ratio for an effective detector area 
of 10 4 m 2 in a year. 



m x (GeV) 


Indirect detection rate (events/10 4 m 2 /year) 


Ns/N B 


50 


2.3 x 10 1 


~ 1 


500 


2 x 10 2 


~ 10 2 


2000 


1.7 x 10 2 


~ 10 4 



In the next three sections we make quantitative calculations of the WIMP annihilation 
rates expected in supersymmetry and supergravity models to replace the simple estimate 
from Eq.(8). 

III. THE MINIMAL SUPERSYMMETRIC MODEL 

A supersymmetry (SUSY) between fermions and bosons is a compelling extension of the 
Standard Model (SM) for several reasons. 

(i) It provides an elegant solution to the fine tuning problem; 

(ii) With the particle mass spectra of the minimal supersymmetric standard model 
(MSSM) the evolution of renormalization group equations is consistent with a grand unified 
scale Mqut ~2x 10 16 GeV and a SUSY mass scale in the range Mz ^ Msusy ~ 1 TeV; 

(iii) With a large top quark Yukawa coupling at a grand unified scale (M G ut), radia- 
tive corrections can drive a Higgs boson mass squared parameter negative, spontaneously 
breaking the electroweak symmetry, and naturally explain the origin of the electroweak scale; 
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(iv) If i?-parityQ is conserved, the lightest neutralino in supersymmetric models can be 
a candidate for cold dark matter; 

The minimal supersymmetric standard model (MSSM) [K0] is the minimal extension of 
the standard model (SM) with (i) a supersymmetric partner for each SM particle and every 
Higgs boson, (ii) two Higgs doublets H\ and H2 such that H\ couples to the fermions with 
£3 = —1/2 and H2 couples to the fermions with £3 = +1/2, (iii) a Higgs mixing parameter 
ft assumed to be real, (iv) three real symmetric 3x3 matrices of squark mass-squared 
parameters, MS, M?. and M~ , (v) two real symmetric 3x3 matrices of slepton mass-squared 
parameters, M~ and M|, (vi) three real 3x3 matrices of trilinear Yukawa parameters, 
Ajj, A D and A E , and (vii) a conserved i?-parity R = (_i)3^+l+25 _ (_^3(b-l)+25_ In 
the MSSM, the lightest supersymmetric particle (LSP) is stable and is a very attractive 
candidate for cosmological dark matter. 

The neutral electroweak gauginos (B and W 3 ) and the Higgsinos (H® and H®) have the 
same quantum numbers. They mix and generate four mass eigenstates called neutralinos. 
The neutralino mass matrix in the basis of (B, W 3 , H®, H®) states is 

— M z cos (3 sin 9 W M z sin (3 sin 9 W \ 
M z cos (3 cos 9 W — M z sin. (3 cos 9 w 
-fi 
-fi 



M 



N 





M 2 



(22) 



^ 2 

—Mz cos (3 sin 9 W Mz cos (3 cos 9 W 
\ Mz sin (3 sin 9 W —Mz sin f3 cos 9 W 

This mass matrix is symmetric and can be diagonalized by a single matrix [OT . The diago- 
nalization is given by the following analytic expressions [p2|,|53| 



/ 



e 1 M x o = -( 1 -a- 1 -C 2 )V> + 



eiM x o = +(l a - \C2) 1/2 



e 3 M Y o 



- 7^) 1/2 " 



3 v 2 6 



6 4 M x o = +(l a - \C2Y' 2 + 



1 !^ 
-a - -C 2 + 



1 l r 



1 !^ 

-a - -C 2 + 



1 1 

-a - -6 2 



(8a 




1/2 








(8a 


- 3O2 


)l/2 




c 3 




(8a 


- 3O2 


)l/2 








(8a 


-§^ 


1/2 



1/2 



+ -(Mi + M 2 ) 



1/2 



+ -(Mi + M 2 ) 



1/2 



+ -(Mi + M 2 ) 



1/2 



+ -(Mx + M 2 ) 



(23a) 



(23b) 



(23c) 



(23d) 



where e« is the sign of the ith eigenvalue of the neutralino mass matrix, and 



3 The R— parity is defined as R = (— i) 3 - B + L + 25 ' i where B = Baryon number, L = Lepton number 
and S= Spin. R = +1 for the SM particles and Higgs bosons, and R = — 1 for their superpartners. 
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C 2 = [M X M 2 -M 2 Z - li 2 ) - -(Mi + M 2 ) 



-(Mi + M 2 ) 3 + J(Mi + M 2 ){M 1 M 2 -M 2 Z - li 2 ) + (Mi + M 2 )li 2 
o z 



+(Mi cos 2 9 W + M 2 sin 2 6 W )M 2 Z - llM 2 z sin 2(3 , 
C 4 = (Mi cos 2 6 W + M 2 sin 2 6 w )M 2 z ll sin 2/5 - M x M 2 ii 2 

+^(Mi + M 2 )[(Mi + M 2 )/i 2 + (M x cos 2 V + M 2 sin 2 W )M| - /iMf sin 

+-1(M 1 M 2 -Ml- M 2 )(Mx + M 2 ) 2 - -^-(Mi + M 2 ) 4 , 
lb 25o 



(24a) 



(24b) 



-Re 



D 



2V3 

_4[/3 _ 27S 2 



S + i(D/27) 1/2 
U = 



1/3 



l -C 2 2 -AC^ S 



2 

"C3 ~ ^C*2 + 



(24c) 
(24d) 
(24b) 



The masses given by the above expressions are usually not in the order M x o < M x o < M x o < 
M x o, but the eigenstates can be relabeled. 

The mixing matrix Z, defined by .Mdiag = ZM. N Z~ l , is given by analytic expressions 

1 M x - €iM x o 



52,53 



Zi 2 

Zn tan W M 2 — ejM x o 



(25a) 



Zi3 li[M 2 - e;M x o] [Mi - e;M x o] - M| sin (3 cos/3[(Mi - M 2 ) cos 2 ^ + M 2 - e;M x o 



^3 

2;1 



Zjl 



and 



M z [M 2 - 6jM x o] sin W [// cos /3 + e;M x o sin /3] 



(25b) 



- ei M x o[M 2 - Ci M o][Mi - Ci M 0] - M| cos 2 /3[(Mi - M 2 ) cos 2 9 W + M 2 - e;M x o 



M z [M 2 - eiM x o] sin W [/i cos + qM^ sin (3] 



l+ ( z _R\ + ( z jA + ( z _* 

\Zi\J \Zi\J \Zii 



1-1/2 



(25c) 



(26) 



In most SUSY models, the lightest neutralino (x?) is the preferred the dark matter 
particle; it is a mixture of gauginos and higgsinos, 



X? = z n B + z 12 W 3 + z l3 H° + z u Hl 
The gaugino fraction of \i can be defined as 



fa = z 2 n + z\ 2 



(27) 



(28) 



and its higgisino fraction as 



Ih = z\ z + z\ A . 



(29) 
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The lightest neutralino is more gaugino-like for /i > M 2 and gauginos annihilate mostly 
into heavy quarks; the lightest neutralino becomes more higgisino-like for // % M 2 and the 
higgsinos annihilate dominantly into gauge bosons. 

The minimal supersymmetric model has 63 free parameters with real mass matrices and 
couplings. To make the analysis of the MSSM tractable, we assume a common value for 
the masses of scalar fermions and the trilinear couplings (Msusy = m* = Af). The most 
relevant of the remaining free parameters are the SU(2) gaugino mass (M 2 ), the Higgs mixing 
parameter the ratio of vacuum expectation values (VEVs) of Higgs fields (tan/3 = v 2 jv\) 
and the CP-odd Higgs-boson mass (m^). 



A. The Effects of Muon Energy Threshold 

The indirect search for SUSY neutralino dark matter with high energy neutrinos is a 
promising approach to detect the annihilations of neutralino dark matter accumulated in the 
Sun |31[| - lf47fl . To obtain realistic rate estimates, we consider the muon detection threshold, 



quark hadronization, and solar absorption. 

For neutrinos from XX ~^ W + W~ — > /U + z/ /t + X the average energy of the muon is 
about m x /4; it is about m x /6 for neutrinos from XX QQ ~ > l^v^O. + X. The detection 
energy threshold for muons makes a significant impact on the indirect detection rate. The 
effect of muon threshold on the indirect detection rate from MSSM neutralino dark matter 
annihilations in the Sun is demonstrated in Figs. 1, 2 and 3 for m x o = 200 GeV, 500 GeV, 
and 1000 GeV, with two values of tan j3 and two sets of M 2 and /i: (i) ji > M 2 such that the 
lightest neutralino is more gaugino-like, and (ii) fi < M 2 such that the Xi is more higgsino- 
like. For simplicity, Mgusy — 1.5m x o is assumed for a common squark mass (rriq), a common 
slepton mass (mi), and a common trilinear coupling (A u = Ad = Ai). We note that the 
higgsino-like neutralino has a substantially higher rate with a higher average muon energy 
than a gaugino-like neutralino. In addition, a higher detector energy threshold effectively 
reduces the indirect detection rate for background and the signal with m x o < 200 GeV. 
For m x o > 500 GeV, the indirect detection rate is not very sensitive to the muon energy 
threshold. 



B. Hadronization of Quarks and Solar Absorption 

There are uncertainties in the signal rate calculation form hadronization of quarks and 
from solar absorption that we now discuss. For a gaugino-like neutralino, the dominant 
annihilation modes are bb or rf for m x < mw± and bb for m x > my/±. For a higgsino-like 



neutralino, the dominant modes are ZZ or W + W~ for m x < m t and ti for m x > m t [f>4 
The effects of hadronization of heavy quarks into heavy mesons or heavy baryons reduces 
the energy of the neutrinos from the semileptonic decays |54 



Moreover, neutrinos produced in the Sun may be absorbed before they escape the solar 
medium. We use the Ritz and Seckel parameterization up to energies around 1 TeV for 
all neutrino producing annihilation modes [BBJ to estimate the neutrino absorption losses. 



Above the TeV scale, extra corrections are treated in a manner similar to work by Edsjo 
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Indirect detectors generally have neutrino energy thresholds of at least a few tens of GeV. 
We have chosen 25 GeV for the muon energy threshold in our analysis which corresponds 
to a neutrino energy of about 50 GeV. 

The combined effects of hadronization, solar absorption and muon threshold energy can 
vary from total obliteration of the neutrino signal for a neutralino mass around 100 GeV to 
merely a reduction by less than one half in neutrino signal for a neutralino mass close to 
1 TeV. Figures 1, 2 and 3 show representative calculations for a typical MSSM neutralino 
dark matter with a mass of 200 GeV, 500 GeV and 1000 GeV, respectively. The effects 
become far less significant for heavier neutralinos. However, higher energy neutrinos have 
a larger solar absorption rate, therefore, slightly more neutrinos from annihilations of 1000 
GeV neutralinos are absorbed than those from 500 GeV neutralinos. 

C. Relic Density and Indirect Detection Rate 

For large tan/5 with heavy SUSY particles, the neutralino annihilation cross section 
is enhanced by the s-channel diagrams involving broad resonance poles of Higgs bosons, 
XiXi ~~ > ^4°; H° — > bb, rf, where A is the CP odd pseudoscalar (A), and H° is the heavier 
CP even scalar (H) . Consequently, the neutralino relic density lies within the cosmologically 
interesting region in much of the parameter space in SUSY models |58|-|63|j . 

In Figure 4, we present relic density of the neutralino dark matter (Q x oh 2 ) versus M 2 
for several values of fi with Msusy taken to be the larger of 300 GeV and 1.5m x o. For 200 
GeV < M 2 < 1200 GeV with 3 < tan/3 < 50 and M SUS y ^ max(300 GeV, 1.5m x o). the relic 
density can be in the cosmological interesting region. 

The indirect detection rate (dNm/dA) in events/km 2 /year is presented versus M 2 in 
Fig. 5. Here we have assumed a detector threshold energy of muons to be 25 GeV. Figure 6 
shows regions in the parameter plane of (M 2 , fx) that may yield an indirect detection rate in 
events/km 2 /year: (i) <iiV ID > 100 (dark shading), (ii) 100 > <tW id > 10 (intermediate) (iii) 
dNiD < 10 (light shading). Parameters in the blank regions do not generate a cosmologically 
interesting relic density (0.05 < f2 x o < 0.3) for the neutralino dark matter. 

The indirect detection rate can be significant (> 10 events/km 2 /year) in large regions of 
the parameter space with 10 ^ tan/3 ^ 50: (i) in the neighborhood of M 2 ~ 500 GeV and 
fi ~ 500 GeV, (ii) in the neighborhood of M 2 ~ 4000 GeV and fi ~ 1200 GeV, (iii) in a 
narrow band with M 2 ^ 1200 GeV and M 2 ~ 2fi, and (iv) in a narrow band with M 2 ~ 400 
GeV and fi > 1200 GeV for tan/3 ~ 50. 

IV. THE MINIMAL SUPERGRAVITY MODEL 

In the minimal supergravity (mSUGRA) model it is assumed that SUSY is broken in a 
hidden sector with SUSY breaking communicated to the observable sector through gravi- 
tational interactions, leading naturally to a common scalar mass (mo), a common gaugino 
mass (777,1/2), a common trilinear coupling (A ), and a bilinear coupling (-Bo), at the GUT 
scale (Mqut ~ 2 x 10 16 GeV). Through minimization of the Higgs potential, the B parame- 
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ter and magnitude of the superpotential Higgs mixing parameter // are related to tan (5 and 
M z- 

The mass matrix of the charginos in the basis of the weak eigenstates (W , H^) has the 
following form 

Mc -{V2M w cosP „ J' (30) 

The sign of the [i term in Eq. ([H]) establishes our sign convention. Analyses of b — > 57 decay 
|6^ , |98|j66H69H69| and muon g-2 measurements [f7D]| - H§9| strongly favor \i > 0. We shall only 



consider [l > and A = 0. 

The SUSY particle masses and couplings at Mz can be predicted by the evolution of 
renormalization group equations (RGEs) from Mqut- The gaugino masses at the weak scale 
(Mz) have the following relations: 

• M l /a 1 = M 2 /a 2 = M 3 /a 3 , 

• Mi ~ 0.4m 1/2 , M 2 ~ 0.8m 1/2 , M 3 ~ 2.7m 1/2 , 

• m^o : m^o : m + : m s ~ 1 : 2 : 2 : 6. 

Al X2 Xi y 

We calculate SUSY particle masses at Mz and the masses and couplings in the Higgs sector 
with one loop corrections in the one-loop effective potential at the scale Q = ^/ m i L m t R 
5?|,|6T]|. At the scale Q, the RGE improved one-loop corrections approximately reproduce 



the dominant two loop corrections p0|-[95|] to the mass of the lighter CP-even scalar. 



We impose the following theoretical requirements on the RGE evolution: 

(i) radiative electroweak symmetry breaking (EWSB) is achieved, 

(ii) the correct vacuum for EWSB is obtained (tachyon free), and 

(iii) the lightest SUSY particle is the lightest neutralino. 

In Figure 7, we present the relic density of the mSUGRA neutralino dark matter (Q x oh 2 ) 
versus mi/2 for several values of mo. For tan (3 ^ 35, cosmo logically interesting relic densities 
are found for m 1( / 2 ^ 500 GeV. For tan/5 < 50, dips in Q x a appear near resonance poles at 
2m x o ~ m H , rriA- 

Figure 8 shows the higgsino fraction of the lightest neutralino in the mSUGRA model. In 
most the of the mi/ 2 , mo plane, Xi is mostly gaugino-like even with mo > 1000 GeV [96,97 



Only in a narrow band when values of mi/ 2 and mo are close to theoretically excluded region, 
can the lightest neutralino have a higgsino fraction larger than 10%. The region excluded 
by the m x + < 103 GeV limit from the chargino search [^] at LEP2 is indicated. 

The indirect detection rate (dNio / dA) in events/km 2 /year versus is presented in 
Fig. 9 for the mSUGRA model. Again we take the detector threshold energy of the muon 
to be 25 GeV. Figure 10 shows regions in the parameter plane of (m^jWo) that may yield 
an indirect detection rate in events/km 2 /year with: (i) 10 > dNn)/dA > 1 (intermediate- 
shading) and (ii) dNio/dA < 1 (light shading). There is no region that has dNm/dA > 10. 
Parameters in the blank regions do not generate a cosmologically interesting relic density 
(0.05 < Q x o < 0.3) for the neutralino dark matter. For tan/3 = 50, the relic density is 
suppressed below the cosmologically interesting value when 2m x o ~ m#, m^. That leads to 
a blank narrow band in Fig. 10(b). 
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V. SUPERGRAVITY MODELS WITH NONUNIVERSAL HIGGS-BOSON 

MASSES AT THE GUT SCALE 



We next consider a supergravity (SUGRA) unified model with non-universal boundary 
conditions for the Higgs bosons at the unified scale (th-gut)- We parameterize the GUT-scale 
Higgs masses asQ 

m Hi (GUT) = (1 + 5i)m = pim , % = 1, 2 (31) 

where 8\ = — 1 and 2 correspond to pi = and 2 The universal mSUGRA model has 
81 = 82 = 0, while 8i = —1 and 8i = 2 correspond to m^(GUT) = and m^(GUT) = 2m . 
The nonuniversality of Higgs-boson masses at mqut can significantly affect the values of 
Higgs masses and couplings at the weak scale. 

In Figures 11 and 12 we present the relic density of neutralino dark matter (Q x oh 2 ) versus 
8\ and 82 for several values of itiq. For tan/3 < 10, the non-universality of the Higgs masses 
makes slight only impact on relic density of the neutralino dark matter. For tan/3 > 35, 
non-universality affects the weak-scale Higgs-bosons masses, the value of the p parameter 
and the neutralino relic density. For tan/3 ~ 35 and 50, dips appear when 2m x o becomes 
close to m H , uia- 

In Fig. 13, we present the higgsino fraction of the lightest neutralino (Jh) as defined 
in Eq. (|29| ) in the plane of (81,82) with mi/ 2 = m = 400 GeV, A = 0. For tan/3 ~ 3, 
an increase in 8\ or 82 generates a larger Higgsino fraction for while for tan/3 ^ 10, a 
decrease in 5\ or an increase in 82 can make fn larger for Xv The higgsino fraction of Xi is 
very sensitive to the value of 82- 

The indirect detection rate (dN^/dA) in events/km 2 /year versus mi/2 is presented in 
Figs. 14 and 15. The threshold detection energy of the muon is assumed be 25 GeV. For 
tan/3 > 35, the indirect detection rate can be large for chosen values of the parameters. 
Figures 16, 17 and 18 show regions in the parameter plane of (m^, m ) that may yield an 
indirect detection rate in events/km 2 /year, (i) <iiV ID > 10 (dark shading), (ii) 10 > <iiV ID > 1 
(intermediate shading) (iii) dN\D < 1 (light shading). Three cases are considered: (a) 
8 1 = -0.5 and 8 2 = (Fig. 16), (b) 5i = and 8 2 = 0.5 (Fig. 17), and (a) 81 = -0.5 
and 82 = 0.5 (Fig. 18). Parameters in the blank regions do not generate a cosmologically 
interesting relic density (0.05 < Q x o < 0.3) for the neutralino dark matter. For tan/3 ~ 50, 
non-universality can significantly enhance the indirect detection rate for the neutralino dark 
matter. 

For tan (3 ~ 10, the predicted indirect detection rate is interesting only in a very small 
band near the theoretically excluded region with m 1( / 2 ^ 600 GeV and m ~ 3m 1 / 2 - For 
tan/3 ~ 50, the indirect detection rate is significant in a large region of parameter space 
with mi/j < 1000 GeV and m < 800 GeV. 



In Refs. I IOC , 101 1 , m^.(GUT) are parameterized to be (1 + 8j)m, 
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VI. PROSPECTS FOR INDIRECT DETECTION 



For a very large volume second generation neutrino detector with superb angular resolu- 
tion and large statistics for atmospheric neutrinos of order 10 4 events per year for IceCube, 
this background can be reduced, in fact, essentially eliminated. Because the atmopheric 
neutrino background will be known with high precision, the number of events in the angular 
bin of the Sun can be determined by interpolation from adjacent bins. It should be possible 
to estimate the rate of atmospheric neutrinos to the level of a few events, and to achieve a 
background much smaller than the fluctuations on the events in the bin. In such a detector, 
the dominant background comes from neutrinos produced by cosmic ray interacting in the 
Sun's corona. This background is between 10 and 5 events per year per square kilometer for 
muon thresholds between 25 and 100 GeV [Q. 

The next generation high energy neutrino telescopes such as IceCube and ANTARES 
with effective area on the order of 1 km 2 may have the capability to make the first detections 
of supersymmetric cold dark matter through the detection of neutrinos produced via anni- 
hilation of heavy neutralino in the Sun. Considering muon thresholds ranging from 25 GeV 
to 100 GeV (reflecting a variety of detectors and their variation in design), our calculations 
show that a significant fraction of the MSSM parameter space with tan/? > 10 should pro- 
duce a muon event rate above backgrounds in such detectors. The prospects for searching 
mSUGRA parameter space are less optimistic than in the MSSM. If non-universal boundary 
conditions are considered, the prospects of detecting SUGRA neutralino dark matter with 
such detectors are increased. 

We have focused on the flux of muons (0^) generated by high-energy neutrinos (in units 
of 'events/km 2 /year'). The event rate for a thin detector such as the AMANDA can be 
obtained from M x A e g, where A e g is the effective area. For detectors of a volume > 1 km 3 
such as the IceCube, the thin detector approximation is valid for muon energies larger than 
a few hundred GeV. For lower energy muons, the event rates will be higher than those 
obtained from thin detector approximation |43l . 



The IceCube detector will have good angular and energy resolution for high energy muons 
[H|. The measurement of muon angular distribution can be used to focus the search window 
around the direction from the Sun and to reduce the background effectively. The atmospheric 
neutrino background in a solar pixel will be determined accurately. The significance of 
the neutralino signal at the IceCube will be greatly enhanced with suitable background 
subtraction. 

The question of naturalness and fine tuning arises for models with heavy neutralinos. 
Recently, it has been shown that TeV scale values of m can be natural [j99fl . The purpose of 
this paper does not lie with this issue and we put forward our results only as a scan of the 
described parameter space setting aside the subjective criterion of naturalness requirements. 
However, because of naturalness, we do not consider the super-heavy models in which the 
lightest neutralino has mass of order several TeV or greater. 



VII. CONCLUSIONS 

High energy muons produced by neutrinos from relic neutralino annihilation in the Sun 
can lead to promising signals in the ice and the water detectors of high energy neutrinos. 



14 



We find that: 



(i) The effects of detector threshold, hadronization and solar absorption can be essential 
in the evaluation of indirect detection rates, especially for neutralinos with mass below 
a few hundred GeV. 

(ii) In large portions of the MSSM parameter space with tan/3 ^ 10, the indirect detection 
rate is predicted to be greater than ten events/km 2 /year. However, the indirect search 
for neutralino annihilation may be difficult for tan (3 < 10 or m x < 200 GeV. 

(iii) The indirect search for the mSUGRA neutralino dark matter will be challenging. 
Only several years of observation with a square kilometer detector will likely lead to 
a discovery. SUGRA models with non-universal boundary conditions and tan (3 > 35 
give more interesting rates. 

(iv) For large values of large tan f3 the neutralino annihilation cross section and the indirect 
detection rate are enhanced by the s-channel diagrams involving Higgs bosons, XiXi ~ * 
A , H° — > bb, tt, where A is the CP odd pseudoscalar, and H° is the heavier CP even 
scalar, and high energy neutrinos are produced in the b and r decays. 

(v) For the AMANDA experiment, with a muon energy threshold of 25 GeV the back- 
ground muon flux from atmospheric neutrinos in a pixel containing the neutralino 
signal is about 4.6 events/(10 4 m 2 )/year. This background flux is larger than the 
signal rate in most regions of the parameter space in SUSY models. To improve the 
significance of signal over background, a higher muon energy threshold will be advan- 
tageous if m x > 500 GeV. 

(vi) For the IceCube and the ANTARES experiments, the large number of background 
events will make background subtraction possible with only a small remaining irre- 



ducible background ||3| and allow measurement of the neutralino annihilation signal. 



(vii) SUSY neutralino dark matter with a mass larger than about 200 GeV offers great 
promise for indirect detection experiments. Together with direct detection experiments 
of neutralino dark matter and accelerator experiments at the upgraded Tevatron and 
the CERN Large Hadron Collider, high-energy neutrino telescopes will be able to 
survey large regions of parameter space beyond present experiments. 
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FIG. 1. Indirect detection rate dNm/dA in events/km 2 /year versus the muon detection energy 
threshold for m x o = 200 GeV and M SU SY = 400 GeV with m A = 175 GeV and (a) tan/3 = 3, 
H = 700 GeV, M 2 = 400 GeV, and y = 250 GeV, M 2 = 485 GeV, as well as (b) tan p = 35, 
y = 700 GeV, M 2 = 400 GeV, and y = 250 GeV, M 2 = 445 GeV. Also shown is the 3a signal rate 
(= 3\/Nb), where Nb is the number of background events from the atmospheric neutrinos. 
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FIG. 2. Indirect detection rate dNm/dA in events/km 2 /year versus the muon detection energy 
threshold for m x o = 500 GeV and Msusy = 750 GeV with niA = 500 GeV and (a) tan/3 = 3, 
H = 1400 GeV, M 2 = 1000 GeV, and fi = 520 GeV, M 2 = 1280 GeV, as well as (b) tan/3 = 35, 
fi = 1400 GeV, M 2 = 1000 GeV, and fx = 520 GeV, M 2 = 1160 GeV. Also shown is the 3a signal 
rate (= 3^/Nb), where Nb is the number of background events from the atmospheric neutrinos. 
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FIG. 3. Indirect detection rate dNm/dA in events/km 2 /year versus the muon detection energy 
threshold for m x o = 1000 GeV and M SUSY = 1500 GeV with m A = 1000 GeV and (a) tan (3 = 3, 
H = 1040 GeV, M 2 = 2060 GeV, and p = 1860 GeV, M 2 = 2000 GeV, as well as (b) tan (3 = 35, 
H = 1040 GeV, M 2 = 2060 GeV, and \i = 1500 GeV, M 2 = 2000 GeV. Also shown is the 3a signal 
rate (= 3V Nb), where Nb is the number of background events from the atmospheric neutrinos. 
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FIG. 4. Relic density of the neutralino dark matter (Q x oh 2 ) versus M% in the MSSM with 
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tan/3 = 35 and (d) tan/3 = 50. 
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FIG. 5. Indirect detection rate of the MSSM neutralino dark matter from the Sun (dNiv/dA) 
in Events/km 2 /year versus M2 with Mstjsy = MAX(300 GeV, 1.5m 0) and several values of // for 
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FIG. 6. Regions of indirect detection rate of the MSSM neutralino dark matter in the (M2, //) 
plane with M SU sy = MAX(300 GeV, 1.5m x o) for (a) tan/3 = 10 and (b) tan/3 = 50. The region 
with dark shading has dNm/dA > 100; the region with light shading has dNm/dA < 10; and the 
region with intermediate shading has 100 > dNm/dA > 10 (event /km 2 /year). The blank regions 
do not have a cosmologically interesting relic density (0.05 < fi^o < 0.3) for the neutralino dark 
matter. 
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FIG. 7. Relic density of the neutralino dark matter (fi^o/i 2 ) versus mi/2 in the mSUGRA 
model with (j, > 0, Aq = and several values of mo for (a) tan/3 = 3, (b) tan/3 = 10, (c) tan/3 = 35 
and (d) tan/3 = 50. 
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FIG. 8. Contours of Higgsino fraction for the neutralino dark matter of minimal supergravity 
in the (m-i/2; m -o) plane for (a) tan/3 = 10 and (b) tan/? = 50. Also shown are the parts of the 
parameter space (i) excluded by theoretical requirements, or (ii) excluded by the chargino search 
at LEP 2. 
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Indirect detection rate of the neutralino dark matter from the Sun (dNn)/dA) in 
/year versus m^j with fj, > 0, Aq = and several values of tuq for (a) tan/3 = 3, (b) 
(c) tan (3 = 35 and (d) tan (3 = 50. 
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500 1000 1500 500 1000 1500 2000 
m 1/2 (GeV) m 1/2 (GeV) 

FIG. 10. Regions of indirect detection rate of the mSUGRA neutralino dark matter in the 
(m 1 / 2 ,"T-o) plane with fi > and Aq = for (a) tan/3 = 10 and (b) tan/3 = 50. The region with 
light shading has dN\Y>/dA < 1; and the region with intermediate shading has 10 > dNm/dA > 1 
(event/km 2 /year). There is no region with dark shading that has cL/Vid jdA > 10. The blank regions 
do not have a cosmologically interesting relic density (0.05 < Q x o < 0.3) for the neutralino dark 
matter. Also shown are the parts of the parameter space (i) excluded by theoretical requirements, 
or (ii) excluded by the chargino search at LEP 2. 
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FIG. 11. Relic density of the neutralino dark matter (O^o/i 2 ) versus b\ in a SUGRA model 
with non-universal Higgs masses at Mqut, with 62 = 0, fi > 0, Aq = 0, and several values of too, 
for TO1/2 = 250 GeV with (a) tan/3 = 3 and (b) tan/3 = 10, as well as for mi/2 = 400 GeV with 
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FIG. 12. Relic density of the neutralino dark matter (O^o/i 2 ) versus 62 in a SUGRA model 
with non-universal Higgs masses at Mqut, with 5\ = 0, fj, > 0, Aq = 0, and several values of too, 
for TO1/2 = 250 GeV with (a) tan/3 = 3 and (b) tan/3 = 10, as well as for mi/2 = 400 GeV with 
(c) tan/3 = 35 and (d) tan/3 = 50. 
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FIG. 13. Contours of the Higgsino fraction of the neutralino dark matter (/#) in the ((Ji,^) 
plan within the framework of a sugravity unified model with non-universal Higgs masses at MquTj 
H > 0, mi/2 = m = 400 GeV, for (a) tan/3 = 3, (b) tan (3 = 10, (c) tan/3 = 35 and (d) tan (3 = 50. 
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FIG. 14. Indirect detection rate of the SUGRA neutralino dark matter from the Sun in 
Events/km 2 /year versus m^2 with /i > 0, Ao = 0, and three choices of non-universal boundary 
conditions for the Higgs masses at Mqut (i) <5i = —0.5 and 82 = 0, (ii) 8\ = +0.5 and 82 = 0, (hi) 
5 1 = -0.5 and 5 2 = +0.5, for m = 200 GeV with (a) tan (3 = 3, (b) tan (3 = 10 and (c) tan (3 = 35, 
as well as for mo = 1000 GeV with (d) tan (3 = 50. Also shown is the indirect detection rate in the 
mSUGRA with Si = 5 2 = 0. 
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FIG. 15. Indirect detection rate of the SUGRA neutralino dark matter from the Sun in 
Events/km 2 /year versus mu2 with /i > 0, Ao = 0, and three choices of non-universal boundary 
conditions for the Higgs masses at Mqut (i) cTi = and 62 = —0.5, (ii) 5\ = and 82 = +0.5, (hi) 
5 1 = -0.5 and 5 2 = +0.5, for m = 200 GeV with (a) tan (3 = 3, (b) tan (3 = 10 and (c) tan (3 = 35, 
as well as for mo = 1000 GeV with (d) tan (3 = 50. Also shown is the indirect detection rate in the 
mSUGRA with S 1 = 5 2 = 0. 
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FIG. 16. Regions of indirect detection rate of the neutralino dark matter in the {m\/2, ffio) 
plane of a non-universal SUGRA model with \x > 0, Aq = and non-universal boundary conditions 
S\ = —0.5 and 62 = 0, for (a) tan/3 = 10 and (b) tan/3 = 50. The region with dark shading has 
dNio/dA > 10; the region with light shading has dNm/dA < 1; and the region with intermediate 
shading has 10 > dN^/dA > 1 (event /km 2 /year). The blank regions do not have a cosmologically 
interesting relic density (0.05 < Q x o < 0.3) for the neutralino dark matter. Also shown are the 
parts of the parameter space (i) excluded by theoretical requirements, or (ii) excluded by the 
chargino search at LEP 2. 
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FIG. 17. Regions of indirect detection rate of the neutralino dark matter in the (mi/2)"Jo) 
plane of a non-universal SUGRA model with \x > 0, Aq = and non-universal boundary conditions 
5± = and 62 = +0.5, for (a) tan/3 = 10 and (b) tan/3 = 50. The region with dark shading has 
dNm/dA > 10; the region with light shading has dNm/dA < 1; and the region with intermediate 
shading has 10 > dNi^/dA > 1 (event /km 2 /year). The blank regions do not have a cosmologically 
interesting relic density (0.05 < Q x o < 0.3) for the neutralino dark matter. Also shown are the 
parts of the parameter space (i) excluded by theoretical requirements, or (ii) excluded by the 
chargino search at LEP 2. 
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FIG. 18. Regions of indirect detection rate of the neutralino dark matter in the {m\/2, ffio) 
plane of a non-universal SUGRA model with \x > 0, Aq = and non-universal boundary conditions 
S\ = —0.5 and 82 = +0.5, for (a) tan/3 = 10 and (b) tan/3 = 50. The region with dark shading has 
dNio/dA > 10; the region with light shading has dNxo/dA < 1; and the region with intermediate 
shading has 10 > dNm/dA > 1 (event /km 2 /year). The blank regions do not have a cosmologically 
interesting relic density (0.05 < O^o < 0.3) for the neutralino dark matter. Also shown are the 
parts of the parameter space (i) excluded by theoretical requirements, or (ii) excluded by the 
chargino search at LEP 2. 
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